OZ-93/22, UM-P-93/89, OITS 519 
CP Violation in Fermion Pair Decays of Neutral Boson Particles 

Xiao-Gang He^ J. P. Ma^ and Bruce McKellar^ 
^Institute of Theoretical Science 

• University of Oregon 

. Eugene, OR97403, USA 

and 

■ '^Research Center for High Energy Physics 

. 

School of Physics 

• University of Melbourne 
^ ■ Parkville, Vic. 3052 Australia 



o 



X 



(October, 1993) 



Abstract 



We study CP violation in fermion pair decays of neutral boson particles with 

spin or 1. We study a new asymmetry to measure CP violation in ry, Kl — > 

H^fi^ decays and discuss the possibility of measuring it experimentally. For 

the spin-1 particles case, we study CP violation in the decays of J/^l^ to SU (3) 

octet baryon pairs. We show that these decays can be used to put stringent 

constraints on the electric dipole moments of A, S and H. 
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I. INTRODUCTION 



CP violation has only been observed in the neutral Kaon system [![]. In order to isolate 
the source (or sources) responsible for CP violation, it is important to find CP violation 
in other systems. In this paper we study CP violation in fermion pair decays of a neutral 
boson particle, which is a CP eigenstate and has spin-O(S) or spin-1 (V). 

a) 5 ^ // 

The most general decay amplitude for S decays into a pair of spin-1/2 particles // can 
be parametrized as 

M{S ff) = Uf{pi, si){as + il-,hs)vf{p2, S2) , (1) 

where as and 65 are in general complex numbers. If both as and bs are nonzero, CP is 
violated. One can define a density matrix R for the process S ff, where the /(/) is 
polarized and the polarization is described by a unit polarization vector Si(2) in the /(/) 
rest frame. With the amplitude in eq.(l) the CP violating part of the density matrix in the 
rest frame of S is given by 

Rcp = Nf{lm{asb*s)p ■ (si - S2) - Re{asb*s)p ■ (si x S2)} , (2) 

where Nf is a normalization constant, and p is the three momentum of the fermion /. Rcp 
contains all information about CP violation in the decay. The expectation values of CP odd 
and CPT odd observables are proportional to Im{asb*g). This CP violating parameter can 
be measured by the asymmetery 

where N+(-) indicates the decay events with Si ■ p > 0(<)0. In terms of the parameters in 
the decay amplitude, 

Im{asb*s) 



A{S) = (3f 

- ■ -Js 



Pj\as\' + \b^^' 
P}MsIm{asb*s) 
SttFj 



2 



where T j is the decay width for S ff, (3f = y'l — Arn?j/Ms, and nif and Ms are the 
masses of the fermion /(/) and the scalar S, respectively. The asymmetry A has been 
studied extensively for ti.Kl yU'*'/^" decays Here we will instead study a CP odd 

and CPT even observable which can provide additional information about CP violation. It 
is related to Re{asb*g). We construct the following asymmetry B to probe this CP violating 
parameter, 

where A^"*" and indicates the decay events with (si x S2) • p > (<)0. In terms of the 
parameters in the decay amplitude, we have 

TT Re{asb*s) 



B{S) = ~-(3f 



4''pj\as\'+\b, 



'S\ 

^ (3}MsRe{asb*s) 

321/ ■ ^ 

To experimentally measure A and B, one must know the polarizations of the fermions in 
the final state. The polarizations can be analysed by using certain decay channels of / and 
/. Assuming that the polariztions of / and / are analysed by the decays / — > fipf) + X 
and / — ^ f'iPf) + X, with density matrices given by 

Pf = l + a/Si ■ p// , 

Pf = ^- «/S2 ■ P/' , (7) 

where are constants, and the hat on the momentum indicates the unit vector in the 
direction of the momentum. Using this information, we define a more convenient asymmetry, 

A(S) = Z± — = afA(S) 

B{S) = = -ajajBiS) , (8) 

where and N~^^~^ indicate events with p/' ■ p > (<)0 and (p// x pj,) • pj > (<)0, 

respectively. 



b) - // 

The most general decay amplitude for this decay can be parametrized as 

M{V ^ //) = £^M/(pi)[7;.(a + 675) + (pi^ - P2^.){c + id^^)]vj{p2) , (9) 

where e'^ is the polarization of V and in its rest frame = (0, e). If CP is conserved, d = 0. 
The constants a, b, c and d are in general complex numbers. 

The density matrix for this decay in the rest frame of V, up to a normalization constant, 
is given by 

Rij = [ma(pi, Si)[7i(a + 675) + {pii - p2i){c + Z(i75)]fA(p2, S2) 

X v^{p2,S2)[-fj{a* + 6*75) + {pij-p2j){c* + zc/*75)]ua(pi,Si)] , (10) 

where i and j label three-vector components. 

The CP violating part of this density matrix is given by 

■0 ^* S^h ^ 2iW_ 

= ilad pj\ (si — Sojj (si — So) ■ ppi 

+ imMv(si X S2)i + %^rr—l\ (^i ■ p(p x 82)* - S2 ■ p(p x Si)^)} 

My + /mj 

+ 2ibd*MvPj{s2iSi ■ p - siiS2 ■ p + i(p x (si - 82))^} (11) 
+ 4:icd*MvPiPj{-{si - S2) • p + i(si x 83) ■ p} , 

where My the mass of V. In general, V is produced with polarization, and the polarization 
depends on how V is produced and is different in different productions. However, we can 
construct a similar asymmetry as for the spin-0 decay case, which is independent of the 
polarization due to rotation invariance, to probe the CP violating parameters 0, Re{da*), 
Re{dc*), 

N+ -N- 



B{V) 



N+ + N 



' 96T 



f 



{2mfRe{da*) + (M^ - Amj)Re{dc*)) 



N+ + N- 

-ajafB{V) . (12) 



Here /3 = Jl - 4m^/M^. 



II. CP VIOLATION m S ^FF 



In this section we study the asymmetry B for 77, — * fi^fi^ decays. CP violating tests 
in these systems have been studied before [^,^. AH of them concentrated on the asymmetry 
A. Here we show that the asymmetry B is also a good quantity to study CP violation. It 
reveals information not contained in A. 

Because rj is a pseudo-scalar, if CP is conserved, = 0. The CP violating contributions 
are expected to be small. We can use the decay width to determine 6^. /m6^ is determined 
from — i> 77 ^ yU^/i" with the two intermediate photons on-shell. Using experimental data 
for ?7 — 77 0, one obtains 

|Im6,| = ^^lni±|[647rr(.^ 27) /m,]^/^ 

= 1.59 X 10"^ . (13) 

This amplitude is close to the experimental amplitude determined from Br^rj ^ = 
(5 ± 1) X 10^^ IQ. The real part of the amplitude is |Re6| ~ 0.7 x 10^^. Using these numbers 
we find 

\B{ri)\ = 2 X lO^lRea^ + 2.3Ima^| . (14) 

Here we have assumed that Re6,j and Im6^ have the same sign. The asymmetry A is 5.8 x 
10^(Rea^ - 0.44Ima^). 

The parameter is model dependent. In many model a,, is very small In lepton- 

quark models, the constraint on is from the neutron electric dipole moment. If one 
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assumes no cancellations among different contributions to the neutron electric dipole mo- 
ment, is constrained to be less than 2 x 10^^. However, if one allows cancellations between 
different contributions, it is possible to have relatively large 0. The asymmetry B can 
reach 10"'^ or even larger. The polarization of the muons from 77 — ^ P-'^fJ'~ can be analysed 
hj ^ ez/gZ/^. In this case = 1/3. The rj factory at SACLAY can reach a sensitivity for 
A and B at the level a few % in the near future It may provide interesting information 
about CP violation. 
Kl fi^fi~. 

Using data from Kl — ^ 27 and Ki — > fJ^^fi^, we obtain, |Im6A'| = 2 x 10^^^, and 
|Re6i<-| = (0.14 ± 0.16) x 10^^^. The contributions to the asymmetries A and B from direct 
CP violation are given by 

\A{Kl)\ = 3.6 X \ReaK - 0.07Imai^| , 
\B{Kl)\ = 2 X lO^^lO.lReoi^ + 1.4ImaK| • (15) 

We have used the central values for Rebx and Imbx- In the above analysis we have assumed 
that Kl is a pure CP egeinstate. This is, however, not the whole story. Kl is not a pure 
CP eigenstate, 

Kl = J: {\K2 > +e\K, >) , (16) 

where CP\K2 >= — \K2 >, CP\Ki >= \Ki > and the mixing parameter e is measured to be 
2.27xl0~^e*'^/'^. The asymmetries A and B are related to Im(62(a2+e&i)) and Re(62(a2+e&i)), 
respectively. Here Oj and bi are the amplitudes for Ki . The parameter bi is not 



zero. Using the values for the real and the imaginary parts of bi determined in Ref. |]T0 
and set 02 = 0, we obtain 

5(iri)|a.=o~ 0.3x10-3. (17) 

The CP odd and CPT odd asymmetry A{KL)\a2 ~ 10^'^. If experiments measure larger 
value for A and -B, there must be new physics due to large a2- In many models the paramter 
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02 is predicted to be very small. However there are models which can produce large 02 0. 
B can be as large as 10~^. A Kaon factory may be able to see CP violation. 

In both 7] f^'^fJ'~ and Kl —>■ fj^^fi^ decays, the experimental sensitivities for the 
asymmetry A is slightly better than for B. If Rea is larger than Ima, the asymmetry A is 
a better quantity to measure. However, if it turns out that Ima > Rea, the asymmetry B is 
the better one. 

The decay Bd{Bs) can be used to study CP violation also. However in the 



standard model, the branching ratios for these decays are very small If the standard 
model prediction is correct, it is very difficult to test CP violation using these decay modes. 
One, of cause, should keep in mind that should this decay be discovered with a branching 
ratio larger than the standard model predictions, there must be new physics and large CP 
violation may be observed. The same comments apply to fi^fi^. The same analysis 

can also be carried out for Higgs particle decays [1^ . 



III. V ^ FF 

The decays V ^ ff provide new tests for CP violation. In a previous paper we have 
studied a particular case, J/tp ^ AA, and shown that this is a good place to look for CP 
violation 0. In this section we will carry out a more detailed analysis by including more 
decay channels. 
J/iP B^B^. 

The branching ratio for J/ip decays into baryon pairs i?8 and B^ of the SU{'i) octet is 
typically 10"^. With enough J/ip decay events, we may obtain useful information about 
CP violation. In Eq.(9), the 6-term is a P violating amplitude and is expected to be signif- 
icantly smaller than the P conserving a- and c- amplitudes. We will therefore neglect the 
contribution to the branching ratio from h. The relative strength of the amplitude a and 
c can be determined by studying angular correlations between the polarization of J/ip and 
the direction of momentum. Due to large experimental uncertainties associated with the 
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constants which determine the angular distribution, a and c can not be rehably determined 
separately at present |^. In our numerical estimates we will consider two cases where the 
decay amplitudes are dominated by 1) the a-term, and 2) the c-term, respectively. Assuming 
a and b are real and using the experimental branching ratios compiled by the particle data 
group, we obtain numerical values for the asymmetry B. The results are given in Table I. 

The CP violating d-teim can receive contributions from different sources, the electric 
dipole moment, the CP violating Z -Bg coupling, etc. In the following we estimate the 
contribution from the electric dipole moment dsg of Bg. Here dsg is defined as 

Ld^pole = t^Bsa.^^^BsF^" , (18) 

where F'^" is the field strength of the electromagnetic field. Exchanging a photon between 
Bs and a c-quark, we have the CP violating C--B8 interaction 

2 

Lc-A = ~^]^^^^^P^ ^ V2)ci^,cB^i^^B^ . (19) 

From this we obtain 

erfss • (20) 

^ J/V 

Here we have used the parametrization, < Q\cq^c\J/ip >= e^gy- The value \gv\ is determined 
to be 1.25 GeV^ from J/ip ^ f^~^fJ'~- There are additional contributions to d, for example, 
exchanging a photon between the final Bg and B^. We have checked several contributions 
of this type and find them to be small if electric dipole moment is the only source of CP 
violation. Using the above information, we can express the asymmetry B in terms of the 
electric dipole moment of the baryons. The asymmetry B can be used to put constraints 
on the electric dipole moment. In Table II we give the asymmetry B in terms of the 
electric dipole moment of Bq. Note that because photons are off-shell, dsg is measured at 
= Mj^^ from the measurement of -B. If we assume that the extrapolation follow the 
same dependence as the magnetic dipole moment of Bs, ds^^q^ = Mj^^) is smaller than 
dBsiff = 0)- However, the q^ dependence of the electric dipole moment may be completely 
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different from the magnetic dipole moment. It is possible that does not change very 
much from = to = Mj^^. 

The polarizations of Bg and Bg in J/ip —>■ BgBs, can be analysed by certain decay 
channels of Bg and Bg. There are many decay channels available to carry out such analysis 
1^. The neutron polarization can be analysed hj n peue- The proton polarization may be 
analysed by rescattering. It may be difficulty to carry out such analysis. The polarization of 
A can be analysed by, for example, A pn^ . This decay mode has a large branching ratio 
(64%) and a large (0.642). The polarization of E can also be analysed. For E~, one can 
use E~ — >■ n7c~ . This decay mode has large branching ratio (99.85%) with as- = —0.068. 
The polarization of E° can be analysed by S° A7. This is the dominant decay channel 
for S° (100%). The polarization of S"*" can be analysed by pn^. The branching ratio 

is 51.6% and has a large value for as+ (-0.98). The polarization of H° can be analysed by 
H*^ Att^. This is the main decay channel (100%) and the parameter a^o = —0.411. The 
polarization of S~ can be analysed by Att^, again this decay mode is the dominant 

one (100%) and has a large vaule for a=- (-0.456). 

The asymmetry B may not be useful in providing upper bounds for the electric dipole 
moment for neutron and proton because their electric dipole moments have been constrained 



to be very small, dn < 1.2 x 10 earn []T3[ and dp < 10 ecm [|14[. However useful informa- 



tion about the electric dipole moments for A, S and S can be extracted. The experimental 
upper bound on d/<^ is 1.5 x lO^^^ecm [0]. There are constraints on the strange quark electric 
dipole moment and colour dipole moment from the neutron electric dipole moment dn , 
which follow if one assumes that the contributions to dn do not cancel against each other 



| T5| . It is possible that cancellations do occur for dn but not dA and the constraints from 
dn do not necessarily lead to strong constraints on d\. Alternative experimental approaches 
to d\, such as that presented here, should therefore be pursued. If d\ indeed has a value 
close to its experimental upper bound, the asymmetry B can be as large as 0(10^^). Using 
A — > p-K^ to analyse the polarization, we can obtain B as large as 10^^. With 10'^ it is 
already possible to obtain some interesting results. This experiment can be performed with 
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the Beijing e"*" e~ machine. If 10^ J /if) can be produced, one can improve the upper bound 
on c/a by an order of magnitude. This can be achieved in future J /if) factories. There is 
not much information about the electric dipole moment of S and S. The observable B can 
thus be used to put upper bound on the electric dipole moments of S and H. With 10^ J/iIj 
decays, the sensitivity for the elelctric dipole moment is typically 10~^''ecm. 

Our analysis can also be used for J/ip ^ Assuming that the d-term is mainly due 
to the electric dipole moment di of the lepton, we have 



B = mi- ■ , (21) 

e 4 1 + 2mi/M ' ^ ^ 

where mi is the lepton mass. For J/i/j ^ fJ'~^fJ'~, we have, B = Ax 10~^((i^/10~^^ ecm) which 
may be too small to be measured experimentally. 
T ^ //. 

In principle the asymmetry B can be used to probe CP violation in T — BsBg. For 
these decays the branching ratios are not measured yet. They are smaller than the branching 
ratios for J/ip B^Bs, e.g. Br{T pp) < 9 x 10^^. In order to reach the same sensitivity 
for CP violating parameters as for J/ip B^B^, more T events are needed. It may not be 
practical to study CP violation using these decay modes. 

It may be possible to observe CP violation in T ^ //. One particular interest- 
ing decay mode is T ^ r+r^. The tauon polarization can be analysed by the decays 
T ^ TTU , 2tiv , Svrz/ , evv and ^vv. It has been shown that these decay channels provide 
reasonable sensitivity for tauon polarization analysis [jl6[. Assuming the electric dipole mo- 



ment of the tauon is the source for CP violation in this decay, B = 7 x I0~^dr/{10~^^ ecm). 
The experimetal upper bound on dr is 1.6 x 10^^^ ecm, so the asymmetry B can be as 
large as 10^^. Values of d^- as large as 10~^^ecm can be obtained in model calculations. 
The leptoquark model is one of them. In this model there is a scalar which can couple 
to leptons and quarks. The couplings of the leptoquark scalar to the third generation are 
weakly constrained |T^. It is possible to generate a large dr by exchanging a leptoquark at 
the one loop level. 
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Similar experiments can be carried out for other systems, for example, p ,(f) fi'^jj,' and 
Z —>■ II [|18| . . In particular the factory may provide useful information about CP violation 
in fi^fi'. 

IV. CONCLUSION 

We studied CP violation in fermion pair decays of spin-0 and spin-1 particles using a CP 
odd and CPT even observable. The asymmetry B studied in this paper provides another 
test for CP violation in t], Kl —>■ p^ pT decays. This asymmetry can reveal new information 
which is not contained in the asymmetry A studied previously in the liturature. For spin- 
1 particle case, we studied CP violation in the decays of J/ip to the SU{?>) octet baryon 
pairs. We showed that these decays can be used to put stringent constraints on the electric 
dipole moments of A, S and S. Using the J/ip events accumulated at the Beijing e^e~ 
collider, one may already obtain interesting information about CP violation. We encourage 
our experimental collegaues to carry out such analysis. 
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TABLES 



TABLE I. The asymmetry B. 



Decay Mode 


a-term Dominates (in unit d{GeV)) 


6-term Dominates (in unit d{GeV) 


nn 


3.5 X 102 


9.0 X 10^ 


PP 


3.2 X 10^ 


8.2 X 10^ 


AA 


3.8 X 10^ 


8.4 X 10^ 




3.7 X 10^ 


7.6 X 102 


[I] 


3.1 X 10^ 


6.1 X 10^ 


TABLE IL The asymmetry B in terms of the electric dipole moment of Sg. 


Decay Mode 


a-term Dominates (in unit 10^^/ecm) 


b-term Dominates (in unit 10^^/ecm) 


nn 


1.38d„ 


3.5dn 


PP 


1.25dp 


3.2dp 


AA 


1.48dA 


3.3dA 


SE 


1.46ds 






1.22(h 


2.4f/H 
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